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ABSTRACT 

The emission from individual X-ray sources in the Chandra Deep Fields and XMM-Newton 
Lockman Hole shows that almost half of the hard X-ray background above 6 keV is unre- 
solved and implies the existence of a missing population of heavily obscured active galactic 
nuclei (AGN). We have stacked the 0.5-8 keV X-ray emission from optical sources in the 
Great Observatories Origins Deep Survey (GOODS; which covers the Chandra Deep Fields) 
to determine whether these galaxies, which are individually undetected in X-rays, are hosting 
the hypothesised missing AGN. In the 0.5-6 keV energy range the stacked-source emission 
corresponds to the remaining 10-20 per cent of the total background - the fraction that has not 
been resolved by Chandra. The spectrum of the stacked emission is consistent with starburst 
activity or weak AGN emission. In the 6-8 keV band, we find that upper limits to the stacked 
X-ray intensity from the GOODS galaxies are consistent with the ~ 40 per cent of the total 
background that remains unresolved, but further selection refinement is required to identify 
the X-ray sources and confirm their contribution. 

Key words: surveys - galaxies: active - galaxies: starburst - X-rays: diffuse background - 
X-rays: galaxies 



1 INTRODUCTION 

Most of the extragalactic X-ray background (XRB) has been 
conclusively shown to be the integrated emission from discrete 
sources, in particular, the accretion light from Active Galactic Nu- 
clei (AGN). The 1-10 keV XRB has a spectral slope of Y = 1.4 
with a 5 -10 per cent spread in measurements of the normalisation 
(see e.g. lRevnivtsev et all2003tlDe Luca & Molen dil2004|) , a large 
amou nt of which is due to field-to-field variations I Barco ns et alj 
2000). Emission from sources resolved by Chandra in the broad 
0.5-2 keV and 2-1 keV bands is able to account for 80- 90 per 
cent of the XRB iMushotzkv et alj l2000l : iGiacconi et all l2002t 
lAlexander et all2003btlBauer et all2004l) . leading to some claims 
that the origin of the background has been solved. Recent analysis 
of very hard (> 5 keV) X-ray data, however, as well as growing 
evidence from infrared and submillimetre studies, indicate that a 
substantial number of hard X-ray emitting AGN may remain to be 
found. 

Although the resolved XRB in the broad 2-10 keV range 
is high, this does not imply that the background is accounted 
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for at lOkeV, since the fraction is dominated by counts in the 
2-6 k eV range to which te lescopes such as Chandra are most sen- 
sitive JWorslev et ail2 005). Narrow energy band source-stacking in 
the XMM-Newton Lockman Hole (XMM-LH) and Chandra Deep 
Fields North (CDF-N) and South (CDF-S) indicate that the re- 
solved fraction of the XRB decreases from ~ 80-90 per cent over 
2-8 keV, to ~ 60 per cent over 6-8 keV, and only ~ 50 per cent 
> 8keV IWorslev et ai]|2004 Eoo3) . These analyses derived the 
resolved fraction by summing, in narrow energy bands, the X-ray 
fluxes from all individually-detected sources. The same source list 
was used for each narrow band considered, regardless of whether 
or not a source is explicitly detected in that band (most, but not 
all, so urces were actually detect ed in the soft < 2keV range; re- 
fer to lWorslev et alj|2004l 12005. and references therein). The un- 
resolved background component has a spectral shape that is con- 
sistent with an unresolved population of highly-obscured AGN 
at redshifts ~ 0.5-1.5 and with absorption column densities of 
~ 10 23 -5 x 10 24 cm . Since the sources have not been individ- 
ually detected in X-rays, their intrinsic, de-absorbed luminosities 
are probably < 5 x 10 43 erg s" 1 (rest-frame 2-10 keV), though 
this depends on important assumptions about their space density 
and absorption column densities. 



© 0000 RAS 



2 M. A. Worsley et al. 



The steeply-rising X-ray source number counts in the 
5-10 kcV band show no evidence for flattenin g towards fainter 
fluxes lHasinger et alJl200ll : lRosati et al Jl2002t iBaldi et al]l2002h . 
Since the average spectral slope of sources is known to be a strong 
functi on of flux - becoming progressivel y harder with dec reasing 
flux lUeda et alJll999t i Fiore et aljf2003t I Alexander eTaP l 2003b: 
IStreblvanska et all2004l) - this also suggests a large, undiscovered 
population of faint, hard X-ray sources. 

Locally, highly obscured AGN outnumber their unobscured 
counterparts by a factor of ~ 4:1. The three nearest luminous 
AGN - Cen A, NGC4945 and Circinus - all have column den- 
sities > 10 23 cm" 2 with NCG4945 and the Circinus galaxy be- 
ing Compton thick (Matt et al. 2000). Population synthesis models 
(e.g. Gilli et al. 2001; Franceschini et al. 2002; Gandhi & Fabian 
12003: lUeda et al JI2003: IComastrill20ol) require large numbers of 
highly obscured AGN in ord er to accou nt for the 30 keV peak in the 
XRB spectrum ISetti & Woltieilll989l). If the ratio o f obscured to 
unobscured AGN remains high (e.g. Gilli et al. 2001 predict 10:1) 
out to high redshift, then a large number of faint, obscured sources 
could remain undetected in deep X-ray observations. At redshifts 
2 > 1 recent, careful studies suggest the obscured:unobscured ra- 
tio is indeed high (e.g. Alexander et al. 2005 find ~ 6:1). 

Highly obscured AGN can show little or no soft X-ray emis- 
sion due to photoelectric absorption, particularly if the cover- 
ing fraction of the obscuring material is high and the fraction of 
reflected light is low. Current X-ray telescopes have been able 
to detect many low-redshi ft obscured AGN, as well as exam- 
ples of type II quasars (e.g. lNorman et alfcOOlfstern et all2002l 
iGandhi et alj2004 . yet sensitivity to faint, hard sources is still lim- 
ited. Highly obscured sources with Seyfert luminosities (e.g. the 
Circinus galaxy) are undetec table even in the CD Fs beyond red- 
shifts of 0.2-1 (see Fig. 6 in lBrandt & HasingedEooH) - the dis- 
tance at which the majority of unobscured AGN are found. Submil- 
limetre and infrared obse rvations (see e.g. lAlexander et alll2003al 
l2005t l Treister et al. 2004), may be the key to identifying the highly 
absorbed sources through the reprocessed AGN emission emerging 
at these wavelengths. IR colour selection using Spitzer data should 
be able to disentangle starb urst emission from the re-ra diated AGN 
emission lAlonso-Herrero et all2004lLacv et alj2004l2005l) . 

The Great Observatories Origins Deep Survey (GOODS) 
toickinson et alJ [20031) is an ongoing campaign of multiwave- 
length deep-field observations comprised of existing and planned 
surveys. The programme is being carried out in two fields, each 
~ 10' x 16' in si ze. The GOODS-N field is contained within the 
~ 2 Ms CDF-N j Alexander etai]l2003bl) and includes the Hub - 
ble Deep Field North (Wil liams et alll996l:lFerguspn et all2000h 
The G OODS-S field lies inside the ~ 1 Ms CDF-S iGiacconi et alJ 
120021) and contains the Hubble Space Telescope (HST) Ultra Deep 
Field observation. As part of the GOODS project the Spitzer Space 
Telescope has also carried out deep infrared surveys in the region 
from 3.6-24 p.m (PI M. Dickinson) and the HST Advanced Camera 
for Surveys (ACS ) has provi ded deep imaging of the fields in four 
broad wavebands iGiavalisco et all2004l) . not to mention a host of 
observations from ground-based telescopes. 

We have performed X-ray stacking of optically-detected 
GOODS galaxies that have not been individually detected in X- 
rays: we specifically ignore all the X-ray sources that have already 
been individually detected in the CDFs - stacking analyses of these 
objects are covered in our earlier work (see Wo rsley et all 2005). 
Since the GOODS survey probes optical galaxies to low luminos- 
ity and high redshift, once the already-detected X-ray sources are 



removed, the remaining 'normal' galaxies are likely candidates to 
host the hypothesised missing population of highly obscured AGN. 



2 ANALYSIS 

2.1 X-ray data and optically-detected sources 

The 1.95 Ms CDF-N covers a total of 447.8 arcmin 2 and is larger 
than the GOODS-N optical survey region, which is focused on 
the CDF-N aim-point. 3cr X-ray source detection sensitivities are 
~ 2.5 x 10~ 17 and ~ 1.4 x 10~ 16 erg cm" 2 s" 1 in the 0.5-2 
and 2-8 kcV bands respectiv ely. We use the main point-source cat- 
alogue (see lAlexander et al]|2003bl) which contains 503 sources. 
The 0.94 Ms CDF-S is also larger than the GOODS-S optical re- 
gion with an area of 391.3 arcmin 2 and detection sensitivities 
of ~ 5.2 x 10~ 17 and ~ 2.8 x 10" 16 erg cm" 2 s" 1 in the 0.5-2 
and 2-8 keV bands respect ively. The main p oint-source catalogue 
contains 326 sources lAlexander et all2003bl) . 

The deep optical imaging of the GOODS regions, using the 
HST Advanced Camera for Surveys (ACS) is now available in the 
four non-overlapping photometric filters: F435 W (-B435), F 606W 
(Vfcoe). F775W (1775) and F850LP (z 850 ) iGiavalisco etaill2004l) . 
The 5a limiting magnitudes in these four bands are 27.9, 28.2, 
27.5 and 27.4 respectively (assuming an aperture of 0.5"). We 
use the rl.lz version of the GOODS ACS catalogue 1 which con- 
tains 32,048 and 29,661 optically-detected sources in the entire 
GOODS-N and GOODS-S regions respectively. 

2.2 X-ray image stacking 

Our st acking pr ocedure is similar to that taken bv [ Lehmer et alJ 
(2005) (also see lBrandtetai]|200lllNandraetalJl2002l) . At each 
optical source position we extract photon counts from the X-ray 
images and exposure times from the X-ray exposure maps. We used 
circular extraction apertures where a fixed 3-pixel radius (1.476") 
was found to give the best signal-to-noise ratio (SNR) compared 
to 2-pixel and 4-pixel apertures. To avoid contamination by X-ray 
detected sources we did not stack any optical sources lying within 
three times th e 90 per cent encircled -energy fraction radius of an 
X-ray source lAlexander et all2 003b). We performed our stacking 
procedure in several narrow energy bands: 0.5-1, 1-2, 2-4, 4—6 
and 6-8 keV. 

The background was estimated using a Monte Carlo approach. 
Each stacking position was randomly shifted in RA and dec. up to 
60" away from its original position but was not allowed to over- 
lap with the stacking aperture of an optical source position, nor lie 
within three times the 90 per cent encircled-energy fraction radius 
of an X-ray detected source. The Monte Carlo procedure was car- 
ried out 10,000 times (except during the stacking of all the sources 
in the catalogue, where only 1000 trials were computationally fea- 
sible due to the large number of sources; see Table[3}, and the mean 
and variance of the background level determined. 

We quantify the significance of the detection of X-ray emis- 
sion from the stacked objects in terms of the signal-to-noise ratio 
(SNR). This is given by 



where S is the total number of stacked-source counts (i.e sources 

1 Available at http://www.stsci.edu/science/goods/ 
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Table 1. The probability, P{ a i se (cr), of falsely recording a given SNR, 
a, in the stacking signal due to background fluctuations alone, when no 
sources are actually present. The probability of a genuine detection is thus 
1 — Pf a i sc (cr). The values correspond to the one-tailed integrated Gaussian 
probabilities i.e. P ialBC (a) = (2tt)- 1 /2 f°° exp(-x 2 /2) dx. 



SNR (a) P faIse 



0.5 


0.309 


1 


0.159 


1.5 


0.0668 


2 


0.0228 


3 


1.35 x 10~ 3 


i 


3.17 x 10" 5 


5 


2.87 x 10" 7 



and background) and B is the total number of Monte Carlo stacked- 
background counts (i.e. background only). For ease of reference, 
TableQgives the probability of a false detection (i.e. the probability 
of recording the excess counts above background by chance alone), 
for various SNRs. 

Total count-rates were derived by dividing the total number of 
stacked counts by the total stacked exposure time. These were then 
conve rted to fluxes us ing counts-to-flux conversion factors from 
lAlexander etaD (2003b). We assumed a Y = 1.4 power-law plus 
Galactic abs orption of jV H = 1.3 x 10 20 and 8.8 x 10 19 c m" 2 for 
the CDF-N lLockmarl20ol) and CDF-S <Stark et alll993) respec- 
tively. The conversion factors are given in Tableland include the 
necessary corrections for Galactic absorption and the low-energy 
absorption seen in Chandra due to molecular contamination of the 
telescope's optical blocking filters. 

The size of the Chandra point spread function (PSF) varies 
strongly with off-axis angle. The radius enclosing 50 per cent of 
the counts from a point-source is only < 0.5" at the aim-point ris- 
ing to > 4" at an off-axis angle of 10' (at 1.5 keV). Since the sky 
density of the optical sources is high enough to result in overlap of 
the PSFs towards the edges of the fields we restricted our analysis 
to the central 4.5' of the X-ray field to avoid confusion between 
sources, and between sources and background. Stacking sources in 
these central regions also enabled the use of a 3-pixel fixed aperture 
size, which encloses the core of the PSF. The PSF remains fairly 
compact, and more importantly, well-measured over these central 
regions and so the stacked counts can be easily corrected for PSF 
effects. We used the data 2 from a circular parameterisation of the 
PSF in order to calculate and correct-for the encircled-energy frac- 
tion at each aperture position and energy band. The average encir- 
cled energy fractions in each of the energy bands for the CDF-N 
and CDF-S are given in Table|2| 

The total stacked-source fluxes were finally converted to in- 
tensity over the solid angle on the sky in which the sources were 
stacked; in this case this is simply the circular central 4.5' region 
(~ 63.6 arcmin 2 ), minus the area of the regions around X-ray de- 
tected sources which we excluded from our analysis. The excluded 
solid angles are 6.6 and 4.3 arcmin 2 in the CDF-N and CDF-S 
respectively. 



2 Enclosed Count Fractions (ECF) using circular apertures based upon 
the SAOsac model of the Chandra PSF (2005 June 24) which is available at 
http://cxc.harvard.edu/cal/Hrma/psf/ECF/hrmaD1996-12-20hrci.ecf.N0002.fits 



2.3 X-ray Background Model 

In order to calculate the fraction of the XRB which can be 
attributed t o the stacked sources we assume the XRB spec- 
trum from IWorslev et alj 120051) . This is a 1-8 keV power- 
law of photon index T = 1.41 and with a lkcV nor- 
malisation of 11.6keVcm ~ 2 s _1 sr _1 keV -1 as observed by 
De Luca & Molendi 1 2004). Above 8 keV the analytical model of 
GruberetalJ <1999t) is used, although renormalised to smoothly 
intercept the 1-8 keV power-law. Below 1 keV we use a steeper 
power-law (r ~ 1.6), which intercepts the low-energy extragalac- 
tic XRB measurements from Ro berts & Warwickl 1200 ll) . 

Whilst the slope of the XRB is well known in the 
1-8 keV range, the normalisation measured by different instru- 
me nts shows 5-10 per cent variations. These m easuremen ts (see 
e.g. IVecchi et alJll999l: ILumb et alJliool iRevnivtsev et aljEobH 
De L uca & Molendil I2004F " are typically obtained over large re- 
gions, typically several square degrees. The variations in XRB nor- 
malisation can be explained in terms of uncertainties in the cross- 
calibration between instruments, and the field-to-field variations 
arising from the discrete nature of the sources making up the back- 
ground. In pencil-beam fields such as the CDFs, which only sam- 
ple < 0.1 cleg 2 , fiel d-to-field variations are at least 10 per cent 
jBarcons etall EoOO). and further variations could be due to true 
'cosmic variance' - real differences in XRB level beyond those 
simply arising due to sampling statistics. F or exam ple, spatial clus- 
tering on the sky (see e.g. lGilli et all20 03. 2005) could be respon- 
sible for additional variations in normalisation, whilst clustering in 
redshift space could create additional variations in XRB spectral 
shape. Error in our XRB spectrum is taken from the measurement 
errors quoted for the difference sources: we have not added any 
additional uncertainty to allow for field-to-field variations. 



3 RESULTS & DISCUSSION 
3.1 Optical source stacking 

Table[3]summarises our stacking analyses of the optically-detected 
GOODS sources. The CDF-N and CDF-S fields were stacked sep- 
arately to provide robustness. Stacking of all the optical sources 
in the fields revealed highly significant (~ 10-30ct) detections in 
the soft 0.5-1 and 1-2 keV bands. Strong results are seen in the 
2-4 keV band (7.3cr and 7.8cr for CDF-N and CDF-S respectively) 
and the broad 2-8 keV range (4.5cr and 6.5<r). Detection of flux 
in the very hard 4-6 keV band is tentative at 2.6a and 2. Oct, al- 
though the combined CDF-N and CDF-S detection is significant 
at the ~ 99 per cent level. No significant detection is found in the 
6-8 keV band with only a 1.1a signal in the CDF-S and a non- 
detection (i.e. SNR < 0) in the CDF-N. We would also note that 
we see no obvious evidence for any spatial extension of the X-ray 
emission in the X-ray images associated with the stacking analysis, 
although a more complete investigation would be needed to pursue 
this further. 

We repeated the stacking analyses using sub-samples of 
the GOODS sources, selecting on the zgso-band, which sam- 
ples the ~ 8300-9500 A wavelength range. The publicly- available 
GOODS catalogue is zgso-band selected. The strongest detection of 
a soft (0.5 — 2keV) X-ray signal occurs for the zg^o < 25 galax- 
ies, whilst the highest signal-to-noise of hard X-ray signals tend 
to occur for the brighter sub-samples, although the detections are 
borderline except for a 4-6 keV signal at 4.3<r for the CDF-N 
with ^850 < 23 sources. There is also hint of a 6-8 keV signal for 



© 0000 RAS, MNRAS 000, 000-000 



4 M. A. Worsley et al. 



Table 2. Counts-to-flux conversion factors for the CDF-N and CDF-S data. The factors convert source count-rate (count s~ x ) to flux (erg cm -2 s — 1 ) and 
were computed assuming a V = 1.4 power-law spectrum inclusive of Galactic absorption. The factors include the Galactic absorption correction and also take 
into account the effects of absorption due to the Chandra molecular contamination. (See Alexander et al. 2003b for details.) The encircled energy fractions 
give the average fraction of source counts lying within the stacking aperture in each of the energy bands for the CDF-N and CDF-S. The values given were 
calculated using for the 'all sources' cases with a 3-pixel (1.476") radius aperture. 



Energy band Counts-to-flux conversion factor Encircled energy fraction 
(keV) (10 — 11 erg cm -2 count -1 ) (3-pixel radius aperture) 

CDF-N CDF-S CDF-N CDF-S 



0.5-1 


0.69186 


0.65327 


0.705 


0.704 


1-2 


0.44369 


0.46358 


0.690 


0.688 


2-4 


1.4583 


1.5830 


0.653 


0.652 


4-6 


2.4763 


2.7689 


0.612 


0.612 


6-8 


8.9967 


10.254 


0.567 


0.567 



Table 3. The results of stacking the emission from optically-detected, but individually X-ray undetected, GOODS sources. For the CDF-N and CDF-S analyses 
the results are shown for each of the different energy bands and also for the different 2850-band magnitude selections applied to the sources included in the 
stacking. If an excess of X-ray counts was detected then the detection SNR is quoted; '-' indicates non-detection, i.e. SNR < 0. 



Energy band SNR (a) 

(keV) 



CDF-N Sources: 




2850 < 21 


2850 < 22 


2850 < 23 


2850 < 24 


2850 < 25 


2850 < 26 


all sources 


Number stacked: 




117 


278 


571 


1160 


2132 


3720 


10,052 




0.5-1 


10.1 


13.8 


15.9 


17.1 


17.1 


16.2 


14.6 




1-2 


10.8 


20.0 


26.4 


28.2 


31.1 


29.8 


27.6 




2-4 


6.8 


8.2 


9.7 


8.4 


7.7 


8.3 


7.3 




4-6 


2.3 


2.9 


1.3 


3.2 


3.1 


3.2 


2.6 




6-8 






1.2 


1.0 


0.4 


0.2 




CDF-S Sources: 




2850 < 21 


2850 < 22 


2850 < 23 


2850 < 24 


2850 < 25 


2850 < 26 


all sources 


Number stacked: 




128 


289 


537 


1034 


1913 


3426 


9599 




0.5-1 


9.4 


9.9 


10.6 


11.5 


11.7 


12.2 


10.1 




1-2 


13.7 


17.8 


20.3 


21.9 


23.4 


23.3 


22.2 




2-4 


3.8 


5.4 


5.6 


7.7 


8.7 


6.7 


7.8 




4-6 


1.7 


1.4 


2.1 


1.7 


2.8 


1.5 


2.0 




6-8 


1.5 


1.9 


0.7 




0.5 


0.7 


1.1 



Table 4. The percentage of the total XRB flux in the signal from stacking the optically-detected, but individually X-ray undetected, GOODS sources. For the 
CDF-N and CDF-S analyses the results are shown for each of the different energy bands and also for the different 2850-band magnitude selections applied 
to the sources included in the stacking. In the case of non-detection (i.e. SNR < 0), a 3cr upper limit is given. The 'already resolved' column indicates the 
fraction of the XRB which has already been accounted for by X-ray detected sources (see lWorslev et 81/2 005 ). 



Energy band Percentage of XRB flux 

(keV) 



CDF-N Sources: 






2 8 50 < 21 


2 850 < 22 


2850 < 23 


2850 < 24 


2 8 50 < 25 


2 8 50 < 26 


all sources 


already 


Number stacked: 






117 


278 


571 


1160 


2132 


3720 


10,052 


resolved 




0.5 


-1 


1.4 ±0.3 


2.8 ±0.4 


4.6 ±0.7 


7.1 ± 1.0 


9.6 ± 1.3 


12.1 ± 1.6 


17.7 ±2.5 


86 ± 11 




1- 


■2 


0.8 ±0.1 


2.1 ±0.2 


4.0 ±0.3 


6.1 ±0.4 


9.1 ± 0.5 


11.5 ±0.7 


17.3 ± 1.0 


83 ± 5 




2- 


4 


1.3 ±0.3 


2.3 ±0.5 


3.9 ± 0.7 


5.0 ±0.9 


6.2± 1.2 


8.8 ± 1.6 


12.6 ± 2.6 


91 ±8 




4- 


-6 


0.9 ±0.6 


1.7±0.9 


3.5 ± 1.2 


3.8 ± 1.8 


5.0 ±2.3 


6.8 ±3.1 


9.2 ± 5.1 


92 ± 13 




6- 


■8 


< 8.9 


< 13.2 


5.2 ± 6.3 


6.4 ±9.2 


3.2 ± 12.4 


1.9 ± 16.4 


< 81.9 


56 ± 12 


CDF-S Sources: 






2850 < 21 


2850 < 22 


2 850 < 23 


2850 < 24 


2850 < 25 


2 8 50 < 26 


all sources 


already 


Number stacked: 






128 


289 


537 


1034 


1913 


3426 


9599 


resolved 




0.5 


-1 


1.9 ±0.4 


3.0 ±0.6 


4.5 ±0.8 


6.6 ± 1.1 


9.0 ± 1.5 


12.7 ±2.0 


17.4 ±3.1 


74 ±9 




1- 


■2 


1.4 ±0.2 


2.8 ±0.3 


4.3 ± 0.4 


6.4 ±0.5 


9.2 ±0.7 


12.1 ±0.9 


19.1 ± 1.4 


72 ± 5 




2- 


■4 


1.1 ±0.5 


2.4 ±0.7 


3.4 ±0.9 


6.4 ± 1.3 


9.7± 1.8 


10.0 ±2.3 


19.4 ±3.8 


83 ± 7 




4- 


-6 


1.0 ±0.9 


1.2 ± 1.3 


2.7 ± 1.8 


3.0 ±2.5 


6.5 ±3.4 


4.7 ±4.4 


10.5 ± 7.4 


82 ± 11 




6- 


■8 


5.4 ±5.1 


9.9 ±7.6 


5.3 ± 10.4 


< 42.1 


6.8 ±18.9 


12.6 ± 25.3 


32.4 ±42.3 


66 ± 11 
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the CDF-N with z 850 < 23 sources (1.2a), and the CDF-S with 
£850 < 22 sources (1.9<r). For the individual sources, we plotted 
X-ray counts against optical flux in various bands but saw no obvi- 
ous correlation, nor was there any obvious trend apparent in similar 
plots of X-ray flux against various optical colours. 

3.2 Source contributions to the unresolved XRB 

Table|4]shows the total flux due to the stacked sources expressed 
as a percentage of the total XRB level (refer to section l2~3l for de- 
tails of the XRB model assumed). Around ~ 15-20 per cent of the 
0.5-4 keV background is contributed by the GOODS sources. In 
the 4-6 keV band the contribution is ~ 10 ± 5 per cent. However, 
only a 3a upper limit of 82 per cent can be placed on the 6-8 keV 
contribution using the CDF-N, with the CDF-S indicating 32 ± 42 
per cent. Both measurements are consistent with both zero and a 
large contribution in this range. 

Fig. Q shows the spectral shape of the stacked X-ray emis- 
sion from the GOODS sources, and the relative contributions to 
the total XRB level. It is encouraging to note that (at least for the 
< 6 keV range) we find the stacked contribution to the X-ray back- 
ground agrees well with the 'missing' intensity once that from X- 
ray detected sources has been considered. The large errors mean 
that there remains some room for other emission, particularly in the 
0.5-2 keV band, even when a ~ 6 per cent contribution from clus- 
ters is taken into account (see e.g. Mo retti et alj2003l) . The residual 
emission can place important constraints on any truly diffuse com- 
ponents lSoltarl2003l) and the total contri bution of accreting black 
holes to reionization jDiikstra et all2004) but it is difficult to place 
any limits on the still-unresolved emission given the large errors 
and without a more thorough understanding of cosmic variance. 

The reduction in sensitivity of Chandra above 6 keV means 
that the instrumental background heavily outweighs any source sig- 
nal and errors in the stacked X-ray intensity become very large. In 
the CDF-N 6-8 keV band the detection of a stacked signal rep- 
resenting ~ 40 per cent of the total XRB - the full unresolved 
intensity in this band - would only correspond to a SNR ~ 2.1a 
(SNR ~ 1.4cr in the CDF-S), with background counts exceeding 
source counts by a factor ~ 100. The 3a upper limit of 82 per cent 
to the CDF-N stacked intensity is easily consistent with what is 
required to account fully for the background. The CDF-S stacked 
intensity gives 32 ± 42 per cent although the error is again very 
large. The 3<r upper limit in the CDF-S exceeds 100 per cent. 

Since GOODS will detect L* galaxies out to redshift z ~ 6 
we would expect the optically-detected sources to include almost 
every galaxy capable of hosting a significant AGN. Although the 
stacked X-ray flux from the GOODS optical sources is consistent 
with the missing XRB intensity, it is still possible that the contribu- 
tion could be much lower and we briefly discuss the possible source 
populations which may be missing from our analysis. If the actual 
contribution is indeed much less than is needed to explain the unre- 
solved XRB, then the missing AGN must come from rare sources, 
occurring with sky densities < lOdeg -2 and so not sampled well 
in the pencil-beam surveys; or, from very faint galaxies that are not 
optically detected in the GOODS survey. 



3.2.1 Rare source populations 

The first of these possibilities requires the missing AGN to oc- 
cur with sky densities < 10 deg -2 , therefore not sampled by the 
~ 0.1 deg 2 deep surveys, and explaining why the stacked source 



emission fails to account for ~ 40 per cent of 6-8 keV XRB emis- 
sion. However, it is also clear that the missing sources must occur 
with sky densities > l deg -2 . This is because to tal XRB emis- 
sion measurements by iDe Luca"& Molendjf 120041) . who analyse 
the background in the fields of ~ 5.5 deg 2 of XMM-Newton ob- 
servations, find the total XRB emission is a T — 1.4 power-law up 
to at least 8 keV. The missing AGN must therefore be present at 
> 1 deg -2 in order to be present in their analysis and provide the 
full XRB flux. 

At densities of 1-10 deg -2 , in order to be responsible for the 
full ~ 40 per cent of unresolved 6-8 keV XRB, the required flux 
per source would need to be > 2 x 10 -13 erg cm -2 s -1 in this 
band. Since only ~ 1 per cent of the 4-6 keV band is unresolved 
iWors lev et alj|20 05). the required sources need to be Compton 
thick with JV H > 2 x 10 24 cm -2 (to have a negligible < 6keV 
contrib ution) and of quasar luminosity (£2-10 kcv ^ 5 x 10 44 ). 
Whilst iGandhi et alJ (2004) have found several type II QSOs 
serendipitously at sky densities of > 10 deg -2 , their sources tend 
to have iV H = 10 23 -10 24 cm" 2 and the contribution to the XRB in 
the 6-8 keV band is only few per cent. 

3.2.2 Optically undetected sources 

The second alternative - that AGN hosts could remain op- 
tically undetected - is supported by the discovery of ex- 
treme X-ray/optical sources (EXOs) which show X-ray fluxes 
of ~ 10 -16 -10 - 15 erg cm -2 s -1 , yet ar e optically undetected 
with z 8 50 > 28 iKoekemoer et alJ 12004 j) . The 7 EXOs in the 
CDF- S have now been de tected using 24 u.m Spitzer observations 
( Koek emoer et all l2004bl) . Interestingly, IWang et alJ 120041) find 
that several of the EXOs show very hard X-ray spectra, consistent 
with type II AGN at redshift < 6 hosted by very under-luminous, 
or very dusty galaxies. It is plausible that X-ray undetected EXOs 
may contribute a few per cent to the total XRB level (those identi- 
fied to date in the CDF-S contribute ~ 0.5-1.5 per cent). 

3.3 Spectral shape: AGN and starburst components 

The softer < 4keV flux seen from the GOODS sources is most 
likely to be dominated by starburst emission in these faint galaxies 
with a photon index T ~ 1.5-2. The average 0.5-2 keV flux of the 
GOODS sources is ~ 2 x 10 -18 erg cm -2 s -1 at a sky density of 
~ 6 x 10 5 deg -2 . A dding this point to a soft X-ray log AT-log S 
diagram (see fig. 4 in |Bauer et all2004l) is consistent with a direct 
extrapolation of the steeply-climbing number-counts distribution of 
star-forming galaxies, which at these very faint fluxes, outnumber 
AGN by an order of magnitude. Starburst activity is a likely ex- 
planation for the < 4 keV stacked emission although weak AGN 
acti vity is also possible. 

IWorslev et alJ 120051) model the missing hard XRB intensity 
as emission from a population of heavily obscured AGN: a simpli- 
fied spectral model consisting of a F — 2 power-law, plus photo- 
electric absorption and a reflection component, is able to explain 
the shape of the missing > 2keV XRB emission inferred from 
the Chandra Deep Fields for redshifts ~ 0.5-1 and column den- 
sities ~ 10 23 - 5 x 10 24 cm -2 (see Fig.gJ. Within the errors, the 
stacked X-ray emission from the GOODS sources does have a spec- 
tral shape which is consistent with that predicted from obscured 
AGN activity. 

Assuming that the missing hard XRB can be explained by the 
hypothesised population of obscured AGN, the difficulty is in re- 
ducing the error in the 6-8 keV stacked Chandra emission to the 
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Figure 1. The extragalactic X-ray background with the resolved and stacked contributions. Panels (a) and (b) shows the data for the CDF-N and CDF-S 
respectively. In each case the solid bl ack line shows the total XRB level and the broken black diamonds show the intensity resolved from individually-detected 
X-ray sources (see Introduction and Worslev et al 2005). The grey diamonds indicate the missing intensity, i.e. the residual remaining once the resolved 
contributions are subtracted from the total XRB level. The solid black crosses now indicate the stacked intensity due to the optically-detected, but individually 
X-ray undetected, GOODS sources determined in the work here. Additionally, for the 4-6 and 6-8 keV bands, 3cr upper limits to the stacked intensity are 
shown. 



point where definite verification or inconsistency can be seen be- 
tween the XRB residual and the stacked GOODS intensity. If, for 
example, only 10 per cent of the GOODS galaxies contain an ob- 
scured AGN, then the overwhelming majority of stacked X-ray po- 
sitions are simply adding background noise to the measurement; 
we need to be able to select the likely AGN candidates. Our at- 
tempt to do this on the basis of the 2350 -band optical magnitude 
was not successful, with no particular 2350 -band magnitude re- 
striction improving the 6-8 keV signal over that obtained when 
all sources were included, although there is a weak improvement 
when only the brighter (£850 % 24) objects are selected. Selec- 
tions of 2850 i$ 23-25 did, however, succeed in increasing the 
SNR of the 4-6 keV band considerably: from 2.6a to 4.3a in the 
CDF-N (z 8 50 < 23); and from from 2.0a to 2.8a in the CDF-N 

(2850 < 25). 

A much better method of discrimination could be the use of 
infrared fluxes which will soon be available from Spitzer observa- 
tions. Highly absorbed AGN are expected to show strong far in- 
frared emission due to the reprocessing of X-ray emission by the 
obscuring dust. The difficulty here would be separating the ob- 
scured AGN candidates from strong starburst galaxies, which are 
also bright in the infrared. 



4 CONCLUSIONS 

We have stacked the emission from X-ray undetected optical 
sources in the GOODS fields to address the question of whether 
they can account for the missing hard XRB implied by the stacking 
of X-ray detected sources l Wor slev et ail 2005). Our central con- 
clusion is 'plausibly, but not certainly, yes'. Whilst X-ray emission 
is detected at high significance in the 0.5-4 keV range, the strong 
downturn in the Chandra sensitivity makes detection more diffi- 
cult at higher energies. By selecting the sources on the basis of 
zsso-band (8300-9500 A) magnitude, we were able to obtain 4.3a 
and 2.8a detections of 4-6 keV emission in the CDF-N and CDF-S 
fields respectively, whilst in the 6-8 keV band we only achieve con- 
straints at the 1.2a and 1.9a levels, respectively. 



The 0.5-4 keV emission represents some 15-20 per cent of 
the XRB with ~ 10 ± 5 per cent in the 4-6 keV band. When added 
to the fractions due to the resolved XRB the totals are consistent 
with 100 per cent, although depending on XRB normalisation (and 
spectral shape < 1 keV), there remains room for contributions from 
other sources. The CDF-S 6-8 keV intensity is 32 ± 42 per cent, 
suggesting that at least half of the unresolved XRB (~ 40 per cent), 
may be lurking in these sources; however, the errors are consider- 
able, and the CDF-N shows no clear excess in the stacked 6-8 keV 
intensity above the instrumental background, although the 3a up- 
per limit to the potential XRB contribution is 82 per cent in the 
CDF-N, which is more than enough to account for the unresolved 
intensity. 

In previous work, we predicted that the missing component 
of the XRB is due to moderate luminosity, highly obscured AGN, 
with absorption column densities of ~ 10 23 - 5 x 10 24 cm -2 , at 
redshift ~ 0.5-1.5. Since GOODS is able to detect L* galaxies out 
to a redshift z ~ 6, we expected it to include the hosts of this heav- 
ily obscured AGN population. The stacked X-ray flux in the CDF-S 
shows a larger contribution in the 6-8 keV band (see Fig.0, pre- 
cisely where the steeply-rising spectrum of these sources becomes 
important. This could be the sought-after signature of the missing 
population, although the poor SNR and the consequently large er- 
rors mean the detection is uncertain. 

Our results are also consistent with a negligible additional 
contribution to the XRB in the 6-8 keV band, which if true, would 
imply that the highly obscured AGN population is hosted by galax- 
ies which are optically-undetected in the GOODS fields. Even with 
the recent discovery of EXOs, it seems unlikely that such a large 
AGN population can be accommodated by such optically faint 
hosts. A strong contribution could also be due to Compton thick 
quasars occurring with sky densities of 1-10 deg -2 . 

In order to boost the signal-to-noise of the stacked GOODS 
signal there remains much work to be done in selecting the appro- 
priate sub-sample of sources, and avoid adding background noise 
to the measurement by including ordinary galaxies. Infrared-bright 
galaxies identified by Spitzer may hold the key to doing this, assum- 
ing differentiation between dusty star-forming galaxies, and poten- 
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tial obscured AGN, can be accomplished. Additional deep field ob- 
servations, particularly further exposure in the CDF-N, could re- 
solve the situation by increasing the sensitivity to individual, faint 
X-ray sources, as well as improve in the signal-to-noise ratio in 
stacked-source signals in the hardest energy bands. Additional in- 
vestigations into EXOs will also provide useful constraints, as will 
shallow wide-area surveys at hard X-ray energies, which can place 
important limits on the contribution from heavily obscured quasars. 

As a final com ment we note that there is goo d agreement with 
the recent work of Hickox & Markevitch (2005) where both the 
resolved sources and unresolved background are measured in the 
0.5-8 keV band using Chandra in the Deep Fields. The total back- 
ground they determine has a normalization 6 per cent less than we 
assume fSection l2~3l . but with an uncertainty of 1 1 per cent. A large 
amount of the difference in normalisations is probably due to dif- 
ferences in analysis method s, the most significa nt of which is our 
'bright end correction' Csee lWorslev et alj20 05) which attempts to 
correct for the poor sampling of the bright end of the X-ray num- 
ber counts distribution. The resolved fractions found by Hickox & 
Markevitch are in good agreement with our results in the bands 
used (1-2 and 2-8 keV). 
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